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Abstract. We study the measurement of the position of atoms as a means to estimate 
the relative phase between two Bose-Einstein condensates. First, we consider N 
atoms released from a double-well trap, forming an interference pattern, and show 
that a simple least-squares fit to the density gives a shot-noise limited sensitivity. 
The shot-noise limit can instead be overcome by using correlation functions of order 
\/N or larger. The measurement of the A^th-order correlation function allows to 
estimate the relative phase at the Heisenberg limit. Phase estimation through the 
measurement of the center-of-mass of the interference pattern can also provide sub- 
shot-noise sensitivity. Finally, we study the effect of the overlap between the two 
clouds on the phase estimation, when Mach-Zehnder interferometry is performed in a 
double-well. We find that a non-zero overlap between the clouds dramatically reduces 
the phase sensitivity. 
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1. Introduction 

Interferometry aims at tlie estimation of the relative phase between two wave-packets. 
In a standard optical interferometer, hke the well known Mach-Zehnder setup [1], these 
two wave-packets correspond to the light traveling inside the two arms of the device, 
and the relative phase 6 is acquired, for instance, as a result of different optical path 
length. After the phase is accumulated, the two wave-packets are recombined through 
a beam-splitter, and the signal at the two output ports depends on 6. The phase can 
be estimated by measuring the difference in intensities between these ports. Apart from 
photons, atoms can be employed for interferometric purposes as well [2j. The atoms 
present some interesting advantages with respect to light, especially due to the non-zero 
mass. In particular, the creation of atomic Bose-Einstein condensates (BECs) opened a 
new chapter in the field of interferometry. The BEC, which behaves like a macroscopic 
matter- wave, constitutes a coherent and well-controllable source of particles. This makes 
the BEC a promising system to measure the electromagnetic [Sj HI [5] or gravitational 
O [71 E] forces. Moreover, the inter-atomic interactions in the BEC are a source of 
nonlinearity, which can be used to create non-classical states [HI [TOl [HI [12] useful to 
overcome the limit imposed by the classical physics on measurement precision [131 [E] • 

A BEC interferometer can be implemented using a double- well trap [151 [HI [13 [HI 
[T9l [201 [211 l22l [23] , where the two wave-packets are localized about the two minima of 
the external potential. In such configuration, a relative phase 6 can be accumulated by 
letting the system evolve in time in presence of an energy difference between the two 
potential minima, and in absence of coupling between the two wells. After this stage, one 
can, for example, recombine the wave-packets by implementing a beam splitter (thereby 
realizing a Mach-Zehnder interferometer). This will imply a further dynamical evolution 
during which atoms oscillate between the wells for a time which must be precisely under 
control, and over which interactions are negligible. 

The recombination of the wave-packets can be done in a simpler way, just by 
releasing them form the double-well trap, so they form an interference pattern, as 
shown in FigiH In this manuscript we discuss how the information about the phase 
can be extracted from this pattern and derive the sensitivity for different estimation 
strategies. The manuscript is organized as follows. In Section [2] we formulate the 
problem and introduce the basic tool - the A^-body correlation function, where A^ is 
the total number of atoms. In Section [3] we demonstrate that by performing a least- 
squares fit to the measured density [15], the estimation sensitivity A9 is bounded by 
the shot-noise. As discussed in detail in Section [H in order to overcome this limit, 
high-order spatial correlation functions must be measured, namely, of order not smaller 
than a/A. In particular, when estimation is performed using the A^-th order correlation 
function, the sensitivity saturates the bound set by the Quantum Fisher Information 
(QFI) [23]. Then, in Section [5] we analyze an estimation scheme based on the detection 
of the position of the center-of-mass of the interference pattern, which can still yield 
sub-shot-noise sensitivity. Finally, in Section [6] we study the sensitivity of the Mach- 
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Zehnder interferometer implemented in a double-well, and we observe that a non-zero 
overlap between the wave-packets dramatically reduces the sensitivity. Some details of 
the calculations are presented in the Appendix. The present manuscript is an extension 
of our previous work [25] . 

2. The model 

To begin the discussion of different estimation methods based on position measurement, 
we introduce the two-mode field operator of a bosonic gas in a double-well potential, 

^(x, t) = ipaix, t)a + ipbix, t)b, 

where creates an atom in the left/right well. With the atoms trapped, the relative 
phase 6 is imprinted between the modes. This stage is represented by a unitary evolution 
tj(9)=e~^^^'^ of the initial state \ipin) of the double-well system. The three operators 

4 = {a}b + Pa )/2, Jy = {a^b - Pa )/2i and = {a^a - Pb )/2 

form a closed algebra of angular momentum. With the phase acquired, the trap is 
switched off and the two clouds described by the mode functions ipa/b{x, t) freely expand. 

The most general quantity, containing the statistical information about the 
positions of particles forming the interference pattern, is the conditional probability 
of finding particles at positions xn = {xi . . .xn)- It can be expressed in terms of the 
A^-th order correlation function Pn{.xn\0) = -^Gn^xn^O), where 

Gn{,Xn, 9) = {^out\^\xi,t) . . . ^\xn, t)^{xN, t)... ^(Xi, t)\^out)- 

Here, \ipout) denotes the state after the interferometric transformation, \ipout) = 
g-«9Jz|^.^^ To provide a compact and useful expression for this probability, we take 
following steps. 

First, we decompose the initial state in the well-population basis, \ipin) = 
Ylin=o^n\n^N—n) and suppose that the expansion coefficients are real and posses the 
symmetry C„ = Cjy^n- As we will argue later, such choice of C„'s is natural in context 
of this work. We switch from the Schrodinger to the Heisenberg representation, where 
the field operator evolves according to, 

^g{x, t) = U^{e)^{x, t)U{e) = ^a{x, t)e'h + t)e-'^b. 

The next setp is to introduce the basis of the coherent phase states [26] defined as 

72^ V J 

(where |0) is the state with zero particles). The action of the field operator on these 
states can be written in a simple form, 

^e{x,t)\^,N) = J^ue{x,ip;t)\ip,N -1), 
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Figure 1. Schematic representation of the interferometric procedure. First, a relative 
phase 9 is imprinted between the wehs. Then, the BECs are released from the trap 
and form an interference pattern. The detectors (symbolically represented as open 
squares) measure the positions of atoms and this data is a starting point for the phase 
estimation. 



where U0{x,ip;t) = ipa{x,t)e2^^^'f'^ + 'ijjh{x,t)e . Next, we expand the Fock states 

in the basis of the coherent states 
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Thus we can easily write the result of action of the field operator on the input state, 
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Now evaluation of G]^{xm,6) (and therefore the pn{xn,6) as well) is straightforward 
- we act times with 4/0 on the input state and calculate the modulus square of the 
result. After normalization we obtain 

2n 2n 



Pn{xn\0) 





TV 

E 



CnCm COS [p (f - n)] COS [ip' (f - m)] 



(2) 



"'-=0 A/(n)(J 

In the remaining part of the manuscript, we fix t large enough so that the interference 
pattern is formed. In this regime, the physical properties of the system change only by 
scaling ~ \/t of the characterisitc dimensions of the system. The probability ([2]) is the 
starting point for the following discussion of various phase estimation strategies. 
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3. Estimation via the fit to the density 

The simplest way of estimating the value of 6 is by fitting the average density to the 
interference pattern, as the position of the maximum depends on the relative phase 
between the two clouds. Such fit is commonly employed with BECs in double-wells, in 
order to determine, for instance, the phase coherence in the system [6l [9l ITS] . 

In the experimental realization, the interference pattern is sampled using M 
bins located at positions Xi {i = 1...M). The number of particles rij in each 
bin is measured m times, giving the set nf'] k = l,...,m. The average occupation 
(ui) = limm_>oo YlT=i nf^"^ /m with free parameter 6 is then fitted to the histogram of the 
measured density {xi,ni},i = 1, ...,M, where fij = Yll^=i'^f^ 1^- ^i^^ e&ch 

of M bins is small, {ui) is given by the average density 

{m) = Gi{x,,e)/\x. (3) 

The value of 9 is determined from the least squares formula 

The fluctuations in each bin, A^n, = limm-s>oo X]fcLi(^!^'' ~" calculated from 

the probability p{ni\0) of detecting Ui particles in the i-th. bin (for details of derrivation, 
see [Appendix A] ), 

p{ni\e) = j dxn.j dxN-n, Pn{xn\0). (5) 

Using Eq.([2]) we obtain 

A^n, = e)^x + [G2(x„ x„ 6) - 6)] (Ax)^. (6) 

In Eq. (j4]), (n,) and A^n^ are assumed to be known, since in the phase estimation stage 
the only measured data are fij. The quantities {ui) and A^n^ are instead constructed 
during the calibration stage, preceding the phase estimation stage, by repeating the 
experiment with different known values of 9. If the number of experiments in the 
calibration is large, and in absence of thermal and technical noise, the measured (nj) 
and A^n, will tend to the theoretical predictions given in Eq. ([3]) and ([6]), respectively. 

Our goal at this point is to determine how the quantum fluctuations A^rij in the 
z-th bin influence the sensitivity of the phase estimation via the fit (jlj). To this end, 
we employ the concept of the maximum likelihood estimation (MLE) [271 [28] . If some 
quantity ^ is measured, the MLE is defined as the choice of 9 which maximizes the 
conditional probability P{^\9) for the occurrence of ^ given 9. That is, the phase 
9 is estimated from the condition -^P{^\9) = 0. In case of the fit discussed here, 
the estimation is based on the measured average occupations rij. If the number of 
measurements m is large, then according to the central limit theorem the probability 
distribution for the average in the i-th bin tends to the Gaussian p{'ni\9) = 

, _ e 2A^ni/m _ every shot the atom counts are correlated between the bins. 

VSvrAni/Y'm 
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Figure 2. Sensitivity of the phase estimation from the fit to the density, as a function 
of \ipin) G A (sohd black hne) plotted using ^/rnA9 with Eg. dTTI) and N = 100 particles. 
The blue dashed line represents the shot-noise limit. The horizontal dotted red line 
indicates the position of the coherent state. Clearly, the sensitivity is bounded by 
the shot-noise. The inset shows the behavior of the sensitivity in the vicinity of the 
coherent state. 



However, in order to link the MLE with the sensitivity of the least squares fit, we 
construct the likelihood function as if the measurement results and nj, with i ^ j, 
were uncorrelated. Thus the total probability of measuring the values {n} = {ni . . . um) 
is a product P{{n}\9) = YlfLiPi^il^)- We note that in this case, indeed, the condition 
for the MLE, ^P{{n}\0) = coincides with Eq.®. The MLE sensitivity saturates the 
Cramer-Rao Lower Bound [271 128], = F~'^. Here F is the Fisher information (FI), 

Therefore, the sensitivity for the least squares fit is given by Eq.Q as well. 

In the following, we demonstrate that this sensitivity is bounded by the shot-noise. 
Let us assume for the moment that the second term in the Eq. ([6]) - which is proportional 
to {AxY - can be neglected. Then, as can be seen from Eq.(l6]), the particle number 
distribution is Poissonian. The FI from ([7]) reads 

A'/ 2 2 

-oo 

with the one-particle probability 

Piim = li\Mx,t)\' + \M^,t)\') + ^{l)Re [ra{x,t)M^,t)e''] . (9) 

We now calculate the Fisher information (|8]) explicitly. As the interference pattern 
is formed after long expansion time, the mode functions can be written as 
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where a = y ip is a Fourier transform of the initial wave-packets, common for ipa 
and ipb and the separation of the wells is 2xo. This gives 

F = mN Hdxl^ (^) 1^^!^, 
j_oo \o-^/ l + acosv9 

with a = j^{Jx) and = 2^§^ + 9. Notice that when the expansion time is long, 
the function varies slowly, as compared to the period of oscillations of siny) and 
cosy). Therefore, in the above expression, one can substitute the oscillatory part with 
its average value. Using the normalization of ip we obtain 



2l + a(-l + v/r^) ^ ' 

As < a < 1, we have < F < mN . Therefore the Fisher information for the fit is 
always smaller than the shot-noise, giving /S.6 > A9sn = '^^/^ two-mode input 

state {A9sN denotes the shot-noise sensitivity). Below we argue that the inclusion of 
the second term in the fluctuations in Eq. ([6]) does not improve the sensitivity. 

In Eq.(l6]), the first term Gi{xi, 6) Ax scales linearly with N while the second term, 
as a function of A^, is a polynomial of the order not higher than two, aN"^ + bN + c. 
The fluctuations A^n, must be positive, thus a > 0. Otherwise, for large N, no matter 
how small Ax, we would have A^rij < 0. If a > 0, the positive second term enlarges the 
fluctuations and thus worsens the sensitivity. When a = 0, the first and second terms in 
Eq.(l6]) scale linearly with iV, and for small Ax the second term can be neglected, thus 
we end up again with Eq.([8])lj. 

In order to calculate the sensitivity A9 = in Eq.( ITT|) . we consider the ground 
states of the BEG in a double well potential. In the two-mode approximation, the 
Hamiltonian of the system reads 

H = -Ejl + ^Jl (12) 

We construct a family of states A by finding ground states of the above Hamiltonian 
for various values of the ratio 7 = 7^^- And so, for 7 > 0, the elements of A are 
number-squeezed states and tend to the twin-Fock state \il)in) = I "f") with 7 — )• 00. 
For 7 < the elements of A are phase-squeezed states [23]. With 7 — —00, the ground 
state of ^ tends to the NOON state \iljin) = -j^i\NO) + \0N)). With 7 = we have a 

coherent state, \tpin) = Notice that for all lipin) G A, the coefficients 

Cn, which were introduced in previous section, are real and symmetric. For each state 
in A, we calculate a = j^{Jx), and insert it into Eq.( !TT]) . The sensitivity shown in Figj2] 
is clearly bounded by the shot noise. 

This limitation for the sensitivity can be explained as follows. The value of the FI 
given by Eq.® is expressed in terms of the single-particle probability. We expect the 
useful non-classical many body correlations to decrease the value of A6', but the FI ([H]) 
is insensitive to these correlations, and thus must be bounded by the shot-noise. In the 

X A similar argument shows that increasing Ax also worsens the sensitivity with respect to Eq.Q 
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next section we demonstrate that the estimation based on the measurement of position 
correlations can improve the phase sensitivity. 

4. Estimation via the correlation functions 

In the estimation protocol discussed in this section, the phase 6 is deduced from the 
measurement of the k-th order correlation function Gk{xk\G)- As previously, we choose 
to deduce 6 using the MLE: a set of k positions is measured, and the phase is chosen 
from the condition -^Gk{xk\9) = 0. After m ^ 1 experiments, = F(ly where 

oo 

N f 1 f d V 

= m-j dx,^-^^ [g^P.ix,\e) J , (13) 

— oo 

with pk{xk\0) = Gk{xk, 0). The coefficient y stands for the number of independent 

drawings of k particles from N , i.e. (^)/ i^Zi)- We notice that by setting = 1, i.e. 
the estimator is a single-particle density, we recover the FI from Eq.(l8]). Therefore, the 
measurement of positions of particles used as independent is, in terms of sensitivity, 
equivalent to fitting the average density to the interference pattern, and is limited by 
the shot-noise. 

Let us calculate the FI for the case k = N, corresponding to the measurement of the 
full A^-body correlation function. We represent the mode functions ipa/b using Eq.lITOl). 
This expression allows to calculate ug{x,ip), and, in turn, the probability (EJ, which is 
then inserted into Eq. (|T3l) . The integrals over space can be performed analytically (see 



Appendix B for details) and the outcome is 

F^,^^)=m-4j2cl[n~-j =m-4A'l = FQ. (14) 

Here, by Fq we denote the QFI, which is a maximal value of the Fisher information with 
respect to all possible measurements [21] . In the case of pure states, the value of the QFI 
is given by 4m times the variance of the phase-shift generator, thus in our case it reads 
Fq = m ■ AA'^Jz- As denoted by open circles in FigJ31 the Eq. flT^ gives A9 = ^/mA9sN 
for the coherent state (7 = 0), and overcomes this bound for all \ipin) G A with 7 < 0. 
The NOON state gives the Heisenberg limit, AOhl — '^j^- 

We now discuss how the sensitivity given by the inverse of Eq.( |T3l) changes for 
k < N. The space integrals for k 1,N cannot be evaluated analytically, thus we 
calculate the FI numerically taking Gaussian wave-packets 

im^m'^y-^ (15) 



0"^/ V vrcr^ 



with the initial width = 0.1 and half of the well separation xq = 1. The Figj3] 
shows how the sensitivity from Eq.( lT3|) for N = 8 atoms changes with increasing k as 
a function of \ipin) G A. The sensitivity improves with growing k, and goes below the 



shot-noise limit at kj^m = 4. 
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Figure 3. The sensitivity y/rnA9 (black solid lines) for iV = 8 calculated by numerical 
integration of the Ea. (fT3|) for various fc as a function of \tpin) S A with 7 < 0. The two 
limits, \/rnlS.9sN and ^fm^Qui^ are denoted by the upper and lower dashed blue lines, 
respectively. The optimal sensitivity, given by the QFI, is drawn with red open circles. 
The inset magnifies the vicinity of the coherent state, showing that the sub-shot-noise 
sensitivity is reached starting from fcmin = 4. 



Therefore, one would have to measure the correlation function of the order of at least 



where N ~ 1000, it would be very difficult to use the correlation function of such high 
order for phase estimation. The biggest difficulty resides in the calibration stage, during 
which one would need to experimentally probe a function of a A; dimensional variable 



In the following section we present a phase estimation scheme based on the 
measurement of the center-of-mass of the interference pattern. Although the probability 
for measuring the center-of-mass at position x is a function of just a one-dimensional 
variable, it can still provide the sub-shot-noise sensitivity. Nevertheless, we will 
demonstrate that the implementation of this estimation protocol can be challenging. 

5. Estimation via the center-of-mass measurement 

5.1. Measurement of all N atoms 

In order to estimate 6 from the measurement of the center-of-mass, one has to go through 
a relatively simple calibration stage. Positions of atoms are recorded independently 
and from this data location of the center-of-mass is deduced. Many repetitions of the 
experiment give the function Pcm{x\0) of a one- dimensional variable. The expression for 
this function can be extracted from the full A^-body probability ([2]) by 



For higher numbers of particles, we numerically checked 





In a realistic experiment with cold atoms. 
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Figure 4. The sensitivity ^/rnA9 (black solid line) for = 100 calculated with 
Eq.p7|) as a function of \ipin) G A with 7 < 0. The values of ^/rnAOsN and ^/rnAOHL 
are denoted by the upper and lower dashed blue lines, respectively. The optimal 
sensitivity, given by the inverse of the QFI, is drawn with the red open circles. 



where "5" is the Dirac delta. To provide an analytical expression for this probability, we 
model the mode-functions by Gaussians as in Eq.f lTS]) . Using a reasonable assumption 
that the initial separation of the wave-packets is much larger than their width, i.e. 
e~^o/'^o <^ 1, we obtain 



Pcm{x\9) 



2a^N 



l + ^(Co + C;vf COS Ud+'^-^x 



^2 



(16) 



The details of this derivation are presented in Appendix C Notice an interesting 



property - the above probability depends on 9 only for states with non-negligible NOON 
components Co and Cat, as already noticed in p9j . 

When Pcm{,x\6) is known, the phase can be estimated using the MLE, as used in 
the previous sections. Then once again the sensitivity is given by the inverse of the FI, 
which can be calculated analytically, 

2 



m 



dx 



Pcm{x\0) \de 



d 



Pcm{x\e) = mN' 



1 



1 



(Co + Cn) 



[17] 



where m is the number of experiments. In FigJUwe plot the sensitivity calculated by the 
inverse of the FI f|T7j) as a function of \ipin) G A with 7 < 0. Although the estimation 
through the center-of-mass is not optimal (A^^ > I — ), the sensitivity can be better 



than the shot-noise, with A^^ — t- A6hl for \ipin) — > NOON. The calibration stage is not 
as difficult as in the case of high-order correlations, however phase estimation based on 
the center-of-mass measurement demands detection of all atoms [301 Ell Ell |33] , as 
we show below. 
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5.2. Measurement of k < N atoms 



If the measurement of the center-of-mass is based on detection oi k < N atoms, the 
probabihty (|T6l) transforms into 

{x\e) 



where Pk{.Xk\0) = J dxN_kPN{xiy\9). The probabihty ffTSl) can be calculated in a manner 



similar to that presented in [Appendix C[ The result is 

{x\e) = 



2a^k 



a cos 



kO H — rr^a; 



a2 



where 



N-k 



i=0 



(iV- A;) 



+k 



(j+fc) ( i ) 



(19) 



(20) 



Notice that for k = N we recover the result from the previous section a = 2CoC 

re using the syi 
f|T9|) can be calculated analytically, 



^{Co + Cn)"^ (as we are using the symmetric states, Cq = Cn)- The FI for the probability 



F = mk'^{l - Vl - a^). (21) 

Let us now evaluate a - and thus F - for various k c:^ N. For k = N and the NOON 
state, we have Cq = Cj\f = giving a = 1 and F = mN"^. From Eq. (l20|) we notice that, 
for any A;, a is the sum oi N — k terms, each depending on the coefficients Ci and Ci+k- 
And so, for k 
and Ci = Cn 
is F 

gives a sensitivity bounded by the shot-noise. Each loss of an atom decreases the FI 
roughly by a factor of N, drastically deteriorating the sensitivity. In Figj5] we plot the 
sensitivity \fm!\Q calculated with the FI from Eq.f l2T|) for various /c ~ iV. To calculate 
a, we choose the subset of l^/'j^) G A which are in the vicinity of the NOON state. The 
Figure shows a dramatic loss of sensitivity as soon as k ^ N . 



N — 1, a will be maximal for a NOON-like state with Cq = C^^i = | 
|. For this state we obtain a = and for large N the value of the FI 
mN. Therefore, the phase estimation using the center-of-mass of — 1 particles 



6. Estimation via the position measurement with the Mach-Zehnder 
Interferometer 



6.1. Formulation 

So far, we focused on the position measurement of atoms released from a double-well 
trap, and studied the phase estimation sensitivity. The fit to the density gives sensitivity 
limited by the shot-noise, and this bound can be overcome by phase estimation with 
correlation functions of the order of at least \/N. As it is difficult to measure these 
correlations in the experiment, it will be challenging to beat the shot-noise limit using 
the interference pattern. Although the sensitivity of the phase estimation based on 
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Figure 5. The sensitivity y^A9 (black solid line) for N — 100 calculated with 
Ea.([2T|) as a function of \ipin) G A with 7 < 0. The values of y^A9sN and y/mAOHL 
are denoted by respectively the upper and the lower dashed blue line. The three 
solid lines correspond to the phase sensitivity for estimation of the center-of-mass with 
different number of particles. For k — 100, the sensitivity is below the shot-noise 
and tends to y/mABm^ for |V'm) — >■ NOON. As soon as /c 7^ iV, the sub-shot-noise 
sensitivity is lost and the value of ^JmAQ increases dramatically. 



the center-of-mass measurement can also be sub-shot-noise, the protocol is extremely 
vulnerable to the loss of particles. 

In the above scenario, the sub-shot-noise sensitivity, which relies on non-classical 
•particle correlations, is reached by directly measuring spatial correlations between the 
atoms forming the interference pattern and using the latter as estimators for the phase 
shift. On the other hand, it is well known that the Mach-Zehdner Interferometer 
(MZI) can easily provide sub-shot-noise sensitivity just by a simple measurement of 
the population imbalance between the two arms and a proper choice of the input state 
I'j/'jn). This is because, in the MZI, the correlations between the two modes carry the 
part of the information contained in the particle correlations which is useful for phase 
estimation. When the clouds are released from the trap and the two modes start to 
overlap, the correlations between the two modes are lost, since an atom detected in the 
overlap region cannot be told to have come from either of the two initially separated 
clouds. This is the reason why it is necessary then to measure directly high-order spatial 
correlation functions in order to reach sub-shot-noise sensitivity. 

It would be thus interesting to quantify the effect of the wave-packets' overlap 
on the sensitivity of the MZI. This analysis has also a practical interest since, in 
the implementation of the atomic MZI, the precision of the population imbalance 
measurement can be improved by opening the trap and letting the clouds expand for 
a while. In this way the density of the clouds drops, facilitating the measurement of 
the number of particles. However, during the expansion, the clouds inevitably start to 
overlap, leading to loss of information about the origin of the particles, as noted above. 
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In this section, we show how the increasing overlap deteriorates the sensitivity of the 
MZI in two estimation scenarios. 

The MZI consists of three stages: two beam-sphtters represented by unitary 
evolution operators e^*^'^"' separated by the phase acquisition e~*^'^^. The atomic MZI 
can be reahzed as follows. Consider a two-mode system governed by the Hamiltonian 
f lT2|) with Ec = 0. The first beam-splitter is done by letting the atoms tunnel between 



the two wells for t 



2 E, 



SO the unitary evolution operator reads Ui = e' 



Then, 



an inter-well barrier is raised, in order to supress the oscillations {Ej = 0) and a phase 
between the wells is imprinted, giving U2 = e"*^"^^. The interferometric sequence is 
closed by another beam-splitter, = e*^"^'". The full evolution operator reads 

U(9) = U3U2U1 = e^^-^-e-^^-^^e"^^-^^ = e-^^■^^ 

where we used commutation relations of the angular momentum operators. 

In order to analyze the sensitivity of the MZI, we introduce the conditional 
probability Piv(xiv|6') of detecting atoms at positions xjy = Xi. . .xn- To evaluate 
this probability for any initial state of the double- well system \ipin)i we take the same 
steps as in Section II. In the Heisenberg representation, the field operator evolves as 

^e{x, t) = U^{e)^{x, t)U{d) (22) 

iPa[x, t) COS \ - ] + Wb[X-, t) Sm 



+ 



V'6(a;,t)cos ( ^ 



1pa{x,t)s\\l [ ^ 



h. 



Again, we express the action of the field operator on \il)in) using the basis of the coherent 
phase-states and obtain Eq.([T]) with 

ue{x,(p;t) = 



ipaix, i) COS ( ^ j + il)b{x, t) sin 



+ 



i)b{x, t) cos ( - j - i)a{x, t) sin ( - 



6 2 



e 2. 



Therefore, the probability Pn{.xn\6) for the MZI is given by Eq.([2]) with the ue{x,ip]t) 
function defined above. 



6.2. Measurement of the population imbalance 

The most common phase-estimation protocol discussed in context of the MZI is the 
measurement of the population imbalance between the two arms of the interferometer. 
In order to assess how the sensitivity of this protocol is influenced by the expansion 
of the wave-packets, we introduce the probability of measuring ul atoms in the left 
sub-space as follows 

00 

PimbinL\0) = (J^^ J dXn^j dXN-nLPN{xN\0). (23) 
-oo 
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This probability depends on the expansion time via ipa,bix, t), which enter the definition 
of pn{xn\6). The sensitivity for various expansion times, if m ^ 1 measurements are 
performed, can be calculated using the error propagation formula, 

1 A^n 



m 



1 ' 

9(n> 



(24) 



where 

= Vimh{nL\Q) 
nL=0 

is the average value of the population imbalance and 



N 



^'^n = Pimb{nL\d) {riL - Yj " (^)' 

nL=0 ^ ^ 



are the associated fiuctuations. The probability fl23|) resembles Eq.([5]), and the above 



moments are calculated as in Appendix A resulting in 

oo oo 

dxGi{x\9)-j and A\ = — - / rffs ^2(^2!^) - 
-00 

The two lowest correlation functions for the MZI read Gi{x\9) = (^^(x, t)) and 

G2{x2\9) = (^'^(xi,t)^J(x2,t)^e(a;2,i!:)^e(a;i,t)), with the field operator from Eq.([22]), 
and the averages calculated with the input state. When the two wave-packets don't 
overlap, i.e. ipa{,x,t)%ljl{x,t) ~ for all a; G M, Eq. flMl) simplifies to 

^ 1 A^4sin^g+(Jf)cos^g 

This is the well-known expression for the sensitivity of the population imbalance between 
separated arms. It gives < A9sn for all I'i/'m) G A with 7 > 0. 

We investigate the impact of the overlap on the sensitivity ( IMl) by modelling the 
free expansion of the wave-packets ipa/b{x,t) by Gaussians, 

and take xq = 1 and the initial width = 0.1. In Fig. |6] we plot the sensitivity 
^/rnA6 taking ^ = and = 100 for three different expansion times r. The initial 
sensitivity deteriorates as soon as the condensates start to overlap, and the sub-shot- 
noise sensitivity is lost for long expansion times. We attribute this decline to the loss 
of information about the correlations between the modes. Therefore, special attention 
has to be paid to avoid the overlap when letting the two trapped condensates spread. 
Although we expect that the expansion facilitates the atom-number measurement, the 
conclusion of this section is that any overlap of the spatial modes has a strong negative 
impact on the sensitivity of the MZI. 
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Figure 6. (a) The sensitivity ^Jmt^Q calculated with Ea. (|24)) for three different 
expansion times r as a function of j'i/'iTi) G -A with 7 > 0. The solid black line 
corresponds to the setup shown in (b), where r = and the wave-packets don't overlap. 
The dashed red line corresponds to (c), where r = 3 and the wave-packets start to 
overlap. The dot-dashed green line corresponds to (d), where r = 10 and the wave- 
packets strongly overlap. Clearly, the sensitivity is influenced by any non-vanishing 
overlap. The values of ^Jml^OsN and ^Jm/S.6HL are denoted by respectively the upper 
and the lower dashed blue line. Here, N ~ 100 and 6 ~ Q. 



6.3. Estimation via the center- of -mass measurement for the MZI 



In Section 15.11 we have demonstrated that when the two wave-packets overlap and form 
an interference pattern, phase estimation based on the center-of-mass measurement can 
give sub-shot-noise sensitivity. Here we study the same estimation strategy apphed in 
the MZI case. Again, we start with the probability Pcm{x\d) of measuring the center-of- 
mass at position x, 



Pcm{x\9) 



dxM S X 



1 \ 
-J^^^ij Pn{xn\0). 
i=l J 



UsingpAr(xAr|^) for the MZI, one can analytically calculate the center-of-mass probability 
only in the limit of small 9, 

N N 

Y,Cffi{x) +eY,CiQ+W{l + 1){N - /) Ui+2{x) - fi{x)) 



Pcm{x\6) 

where 



1=0 



1=0 



(26) 



fi{x) = exp 



ix-Xo{%-l)y 



2al 
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With this probabihty, we can again calculate the sensitivity using the error propagation 
formula [271 128], 

1 A^x 



m 



d{x) 



2 ' 



where 



) oo 

dx pcm{x\9)x and A^x = j dx 



Pcm{x\9)x^ - {xf. 



These two moments can be easily calculated with Eq. fl2^ . giving, in the limit — )■ 0, 



6*^0 



1 

m 



N 



, , (27) 
(J.)2 \x,J 4(J.)2_ 

Notice that when the initial size of the Gaussians tends to zero, we recover the sensitivity 
from Eq.f l25p (in the limit of ^ — > 0). This is not surprising, as when the mode-function 
are point-like, the measurements of the center-of-mass and the measurement of the 
population imbalance are equivalent, and related by Xcm = 2xo;^. Therefore, for small a, 
the measurement of the center-of-mass yields sub-shot-noise sensitivity for all {ipin) ^ 
with 7 > 0. However, for non-zero a, the second term in Eq.( |271) spoils the sensitivity. 



This is because 



N 



> 



is always satisfied. Even if the first term scales better than 



at the shot-noise limit, the other one does not, and will dominate for large N. 

iFiom what we presented in this Section, we conclude that both the population 
imbalance and the center-of-mass measurements can give sub-shot-noise sensitivity for 
the MZI, but both are very sensitive to the growing size of the wave-packets. 



7. Conclusions 



In this manuscript we have discussed in detail how the measurement of positions of atoms 
forming an interference pattern can be useful in context of atom interferometry. We 
showed that the phase estimation based on the fit to the density gives sensitivity limited 
by the shot-noise, because the FI is expressed in terms of the single-particle probability 
only. The sensitivity can be improved below the shot-noise limit by estimating the 
phase using correlation functions of order at least Moreover, we demonstrated that 
the information contained in the A^-th order correlation function allows to perform an 
optimal detection strategy, reaching Heisenberg-limited sensitivity when NOON states 
are used. We also showed that the measurement of the position of the center-of-mass of 
the interference pattern gives sub-shot-noise sensitivity for all states with non-negligible 
NOON component. Both the measurement of high-order correlations and the center- 
of-mass position are difficult to perform. The former requires the construction of a 
function of highly- dimensional variables, while the latter works well only if all atoms 
forming the interference pattern are detected. We attribute the difficulty to obtain the 
sub-shot-noise sensitivity to the fact that, after formation of the interference pattern, 



Phase Estimation from Atom Position Measurements 



17 



the modes cannot be distinguished, and the information useful for interferometry is only 
contained in the correlations between the particles. These correlations are very difficult 
to extract from the experimental data, therefore reaching sub-shot-noise sensitivity with 
two interfering BECs might prove very challenging. In the final part of this work, we 
turned our attention to the MZI, which is known to provide sub-shot-noise sensitivity 
for the simpler measurement of the population imbalance between the two clouds. We 
have shown that the sensitivity of the MZI is strongly infiuenced by a non-zero overlap 
between the two wave-packets, both in case of the population imbalance and the center- 
of-mass measurement. 



Appendix A. Bin fluctuations 

To derive the expression for the average number and the fiuctuations of the atom count 
in a bin, we use the probability from Eq.([5]). With help of Eq.(|2]) we obtain 

r.(r, im / / d^d^' (^\,n,hN-n, C'n^'rn COS [ifi (f - n)] COS (f - m)] 

where 



ag = J dxu*0{x,ip;t)ug{x,ip';t) and bg = J dxu*0{x,ip;t)ug{x,ip';t). 

Axi R-Axi 

Since the distribution of inside the integrals is binomial, one can easily calculate the 
average value as follows 

ni=0 



X 



27r 27r 

CnCm COS [(f (f - ^)] COS (f " "^)] 





N 



n,m=0 



) ( ) 

^ n ) \mJ 

Notice, that using the definition of ag and bg, this can be rewritten as 



{ui) = N JdxiJ dxN-iPN{xN\0) = J dxiGi{xi\6), 

Axi —CO Axi 

where G'i(2;i|^^) = {'^l{xi,t)'$g{xi,t)) is a one-particle density. In the same manner, 
one can demonstrate that the fiuctuations of the number of particles in the i-th bin are 
given by 

A^m = {n'^i) - (ui)^ = j dxi Gi{xi,e) + j j dx2 G 2{x2, 9) - {ui)"^ . 

Axi Axi Axi 

When the bin size is much smaller than the characteristic variation of the mode 
functions, the above integrals are approximated by 

{ui) = Gi{xi,e)Ax and A'^m = Gi{xi,e)Ax + [G2{xi,Xi,e) - Gl{xi,e)] {Axf. 
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Appendix B. The FI for A; = iV 

In this Appendix we prove the Eg. (1141) . We start from the expression for the N- 
particle probability (|2]) for the wave-packets after long expansion time, as in Eq. ffTOj) . A 
straightforward evaluation shows that in this case, the probability can be written as 

PN{SN\9) = [f{xNm\ (B.l) 
where / is real and reads 

n^^m ^ 2f'i t ""-l^-"'^ n m (f^ - M) ■ (-) 

Moreover, we assume that the initial separation of the wells of the trapping potential 
is much larger than the width of the mode functions. Under this assumption, the 
probability fIB.ip is normalized. This probability is inserted into the definition of the 



Fisher information. As / is real, we immediately obtain 



F^k=N) = m-A J dxN[fixN\e)] . 

Now, the order of the integration can be reversed, and the space integrals are performed 
first. As the mode-functions are normalized, the result is following 

F^k=N) = mN2^ E 71^ / / [^(^ - 2^)] T - - I) (cos(^))^-^ 

n=0 \ri> ''^ ^ 

The phase integral can be now easily evaluated, giving 

N / T.T\2 



( N\ 



ra=0 

which is the Eq. ffT^ . 

Appendix C. Evaluation of the center-of-mass probability 

In this Appendix we derive the expression for the probability of detecting the center-of- 
mass at position x, as in Eq. flT6|) . The definition of Pcm{x\9) relates it to the full A-body 
probability by 



Pcm{x\9) 



j c/f iv 5 - ^ ^ pn{xn\0). 



To calculate this probability we assume a long expansion time and use Eqs fIB.ip and 



f lB.2p . Then, we notice that the Dirac delta can be represented as the Fourier transform, 

oo 

Pcm{x\e) = i- jdkjdxi,e-'^^^-^^'^^^^^pN{xM- 



Phase Estimation from Atom Position Measurements 
This equation can be rewritten in a following way 



1 



19 



(c.i; 



where 



Pcm{k\0) 



is the Fourier transform of the probability pcm{x\9). To provide an analytical expression 
for this probability, we assume that the initial wave-packets are Gaussian as in Eq. (fT5|) . 
The integration over space is performed giving 

"^"^ dip [^^ dip' 



Pcmik\0) 
N 



277 Jo 27r 



[I{k,ip,ip')] 



N 



(C.2) 



X 



ECnCm 
COS 



ra,m=0 



\n / \mj 



N 

ip \ n 



cos 



N 



— m 



where 



I{k,ip,ip' 



+ 2 cos 



ip + ip' 



e 2vj2 



with fco = ^ and 



w 



2Nao 



The function I{k,ip,ip') consists of three peaks, located 



at k = ±fco and k = 0. When the well separation 2xo is large compared to the initial 
width of the trapped wave-packets a, then these three peaks are separated and do not 
overlap. Therefore, [I{k, ip, ip')]^ can be approximated by the sum of A^-th powers of its 
three components 



[l[k,ip,ip)\ ~e V 2 Je ^ +e v ^ /e 2^ + 2 cos 



ip -\- ip'\ ]Vfc2 

e 



This result is put into Eq.( ]C.2l) . the phase integrals are evaluated and the result is 



e 2^ + 



Pcm{k\0) 

Then, using Eq. OC.ip . we obtain 



Pcm{x\9) 



we 2JV 



V2TrN 



2w2 _j_ g 



(Co + / Ar/i , ; A 

1 H cos{NO + kox) 



{Co + C^f 



which, with help of definitions of w and /cq, gives Eq.([T6 
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